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Abstract 

We investigate the analytical solution of a new exactly solvable non-central potential of V(r, 9) = 
D (^t^) 2 + r 2 s f n a g + ~i ^fio type, which may be called as the modified non-central Kratzer poten- 
tial. The energy eigenvalues as well as the corresponding eigenfunctions are calculated for various 
values of n and m quantum numbers within the framework of the Nikiforov-Uvarov and Asymtotic 
Iteration Methods for the CO diatomic molecule as an application of this potential. In this paper, 
we first present the effect of the non-central term on the bound-state energy eigenvalues: this effect 
is determined explicitly for different n and m quantum numbers with j3 = 7=0.0, 0.1, 1.0 and 5.0 
values and the results are compared with the findings of the modified Kratzer potential for different 
n and I quantum numbers. Then, we show that the angle-dependent non-central part behaves like 
a centrifugal barrier and it reduces the depth of the attractive potential pocket, which effects the 
bound-state energy eigenvalues. 
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I. INTRODUCTION 



The analytical solution of the radial Schrodinger equation is of high importance in non- 
relativistic quantum mechanics since the wave function contains all the necessary information 
to describe a quantum system fully. There are only a few potentials for which the radial 



Schrodinger equation can be solved explicitly for all n and a 



as super-symm etry (SUSY) 
Method k 



|3j, Nikiforov-Uvarov 



, [111, H2J, H3| and the Pekeris approximation 



11. So far, many methods such 



a, 



14 



8fl, Asymptotic Iteration 



151 ] have been developed and 



applied to solve the radial Schrodinger equation exactly or quasi-exactly for I ^ within a 
given potential. 



The Kratzer or modified Kratzer type potentials [16] we consider in this paper have 
played an important role in the history of the molecular and quantum chemistry and they 
have been so far extensively used to describe the molecular structure and interactions 17]. 
Although such central potentials have been extensively used to describe the properties of the 
system qualitatively, it is known that the dynamical properties of the system should also be 
taken into account for a better description of the system. That means, the potential should 
have not only the radial but also the angle dependent parts. These dynamical properties of 
the system can be taken into account by using a non-central potential and therefore, in this 
paper, we aim to present the effect of the non-central term to calculate the non-zero angular 
momentum solutions of the Schrodinger equation. For this, we take the non-central modified 



Kratzer (is| potential (V(r,6) = D {^f + 



+ 



7 cos ( 



and as an application, we 



calculate the energy eigenvalues of the CO diatomic molecule for various n and m quantum 
numbers. 

The article is organized as follows: In the following two sections, the basic equations 
of the Nikiforov-Uvarov method (NU) as well as the Asymptotic Iteration Method (AIM) 
used so far to solve the resulting differential equation are given briefly. The solutions of the 
Schrodinger equation with the non-central modified Kratzer potential are outlined in the 
fourth section. Then, the solutions of the Schrodinger equation with the non-central modified 
Kratzer potential are obtained by using both methods. The energy eigenvalues and the 
corresponding eigenfunctions are calculated for various values of n and m quantum numbers 
for the CO diatomic molecule and the effect of the non-central term is determined explicitly 
by comparing the results with the findings of the modified Kratzer potential with different 



2 



n and / quantum numbers. The summary and conclusion are provided in section IVIII 



II. NIKIFOROV-UVAROV METHOD 

The NU method is based on the solutions of general second order linear equations with 
special functions. It has been extensively used to solve the non-relativistic Schrodinger 
equation or similar time-independent second-order differential equations and there is an 
extensive literature to look at . However, in order to keep the completeness of 

the paper, we briefly outline NU method here. The idea in the NU method is to convert the 
Schrodinger equation or similar differential equations into the following form: 

^ s y + f M ns) + ^ is) = (i) 

where a(s) and a(s) are polynomials, at most second- degree, and f(s) is a first-degree 



polynomial 



Hence, from Eq. ([I]), the Schrodinger equation or the Schrodinger- 
like equations can be solved analytically by this method. In order to find a particular solution 
of Eq. (TjQ), the following transformation is used: 

iP(s) = <f>(s)y(s) (2) 

it reduces Eq. (pQ) to an equation of hypergeometric type, 

a(s)y" + r(s)y' + Xy = (3) 

and <j)(s) is defined as a logarithmic derivative in the following form and its solutions can be 
obtained from 

(f>'(s)/(f>(s) = n(s)/a{s) (4) 

The other part y(s) is the hypergeometric type function whose polynomial solutions are 
given by the Rodrigues relation 

where B n is the normalization constant and the weight function p(s) must satisfy the con- 
dition 

(ap)' = rp (6) 
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The function ir and the parameter A required for this method are defined as follows 



o — T , a — T 



^(s) = __ ± W (^__ j -~ a + ka , (7) 

A = k + tt' (8) 

On the other hand, in order to find the value of k, the expression under the square root 
must be the square of a polynomial. Thus, a new eigenvalue equation for the Schrodinger 
equation becomes 

A = A n = -nr a (9) 

where 

T (s) = f(s) + 27T(» ( 10 ) 

and its derivative is negative. By comparison of Eqs. (jSJ) and (J9]), we obtain the energy 
eigenvalues. 

III. THE ASYMPTOTIC ITERATION METHOD (AIM) 

In this section we briefly outline the asymptotic iteration method, the details can be 

nnnnn 

found in references [9|, HO, HH, H2J, H3[- The asymptotic iteration method was proposed to 
solve second-order differential equations of the form 

y" = X (x)y' + s (x)y (11) 

where So(x), Ao(x) are functions in Coo(a,b). The variables, so(x) and Ao(x), are sufficiently 
differentiable. 

y" , = X 1 (x)y' + s 1 (x)y (12) 
The second derivative of Eq. ( fTTl) is obtained as 

y"" = X 2 (x)y' + s 2 (x)y (13) 

where 

Xi(x) = X' (x) + s (x) + Xl(x) si(x) = s Q (x) + s (x)X (x) 
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X 2 (x) = \[{x) + si(x) + A (a;)Ai s 2 (x) = s[(x) + s (x)X 1 (x) 
Therefore, for (n + l) th and (n + 2) th derivatives, n—1,2,...., one can get, 

y (n+1) = \ n -i(x)y' + s n ^(x)y 



and 



respectively, where 



t n +V = \ n (x)y' + s n {x)y 



(14) 



(15) 



(16) 



A n (a;) = \n-i_ix) + s n -t(x) + \ (x)\ n -i(x) s n (x) = s n _ x (x) + s Q (x)\ n -i(x) (17) 



The ratio of the (n + 2) and (n + 1) derivatives can be expressed as 



(18) 



For sufficiently large n, we can now introduce the asymptotic aspect of the method; that is, 



A n A n _i 



a 



(19) 



Thus Eq.f fTBl can be reduced to 



(20) 



which yields the general solution of Eq.lflT 



y(x) = exp I — / adt 



C 2 + Ci / exp / A (r) + 2a(r)dr dt 



(21) 



For a given potential such as non-central modified Kratzer, Schrodinger equations are con- 
verted to the form of Eq.tflT]). Then, so(x) and Ao(x) are determined and s n (x) and A n (x) 
parameters are calculated. The energy eigenvalues are obtained by the termination condition 
given by Eq.(Tl9l. 

In this study, we investigate the exact solutions of the Schrodinger equations which the 
relevant second order homogenous linear differential equation takes the following general 
form [lOj, 



ax 



N+l 



[m + 1 



iv, x 



N 



If this equation is compared to Eq. ljTTI) . it entails the following expressions 

, . . / ax N+1 (m + l)\ , . wx N , , 

a, 6 and m are real numbers and w™(N) can be determined from condition Eq. ffT^l) as follows 

m™(-l) = n (2a + 26m + (n + 1)6) (24) 

w™(0) = 2n(2a + 26m + (2n + l)6) (25) 

w™(l) = 3n (2a + 26m + (3n + 1)6) (26) 

w™(2) = 4n (2a + 26m + (4n + 1)6) (27) 

<*(3) = 5n(2a + 26m + (5n + l)6) (28) 
. . . etc 

Hence, these formulae are easily generalized as; 

M fW=KJV + 2)' B ( B+ ^+;> 2 t + 2a ) (29) 

where n = 0,1,2,3,... and N = —1,0,1,2,3,.... The exact eigenfunctions can be derived 
from the following generator: 

y n (x) = C 2 exp (- J a k dt^J (30) 

where n = 0, 1, 2, ... and k > n is the iteration step number. Using termination condition of 
the method given by Eq. (|19|) and Ao and Sq determined by Eq. (|23[) . the eigenfunctions are 
obtained as follows; 

yo{x) = l 

yi{x) = _c 2 (iv + 2v(i-^±^^ +2 

y 2 ( X ) = C 2 (N + 2) V (a + 1) ( 1 - g^±jl + * iP + 2 ) ^ 3 ) x W + » 



a (a + 1) {a + 2) / 36 (p + 3) 



J/3 0*0 = -C 2 — — — 3 1 x 



(N + 2 



„iV+2 



a 



36 2 (p + 3) (p + 4) 2(jV+2 ) 6 3 (P + 3) (p + 4) (p + 5) 3(N+2) 
a(a + l) ' p(p + l)(p + 2) 

. . . etc 

Finally, the following general formula for the exact solutions y n (x) is found as; 

y n {x) = (-1)" C 2 (N + 2) n (a) n 2 F 1 (-n, p + n;a; bx N+2 ) (31) 
wnere [a) n - r((r) , cr - ^ ancl P — (v+2)6 ■ 



IV. ENERGY EIGENVALUES USING NU METHOD 



In this section, we show how to solve the Schrodinger equation for a particle in the 
presence of non-central modified Kratzer potential by using NU method. The standard 
Kratzer potential is defined by V(r) =-D - 2jV Similar to Q, Q, Q, the modified 
Kratzer potential is obtained by adding a D term to the standard Kratzer potential that 
is V(r) = D {j-^) 2 ■ The new exactly solvable non-central modified Kratzer potential we 
examine in this paper is defined as follows: 

VW^dI^Y+V+^L (32) 
\ r / r z sin 9 r z sin 9 

where D is the dissociation energy and a is the equilibrium internuclear separation and (3 
and 7 are strictly positive constants 0, Jj], 18]. The first term of this potential is the modi- 



fied Kratzer potential, the second and third terms are the angle dependent parts. Thus, the 
non-central modified Kratzer potential is defined as a new potential. For the CO diatomic 
molecule, a comparison of the modified Kratzer and the non-central modified Kratzer po- 
tentials is shown in Figure [1] for different f3 and 7 values with 9 = 30° for the non-central 
part. The spectroscopic parameters of the CO diatomic molecule 19] are given in Table [B 
In the spherical coordinates, the Schrodinger equation with the non-central modified 
Kratzer potential is 



2/i 



]_d_ 



,d_ 

Or 



+ 



+ 







r 2 sin 6 89 \ S 09 



d 



+ 



d 2 



D 



r — a 



+ 







+ 



r 2 sin 2 9 d(p 2 
7cos# 



r 2 sin 2 9 r 2 sin 2 9 



ip = Eip 



(33) 



If the spherical total wavefunction as ip (r, 9, <p)— ^^-Y(9, (p) = U(r)Y(9, <p) is inserted into 
Eq. (133|) . the wave equation for the non-central modified Kratzer potential is separated into 
variables and the following equations are obtained: 



d 2 R 2/i 
dr 2 h 2 



E-D 



r-Ar=o 



(34) 



d 2 Q(9) d®{9) 
' + cot 9— ^ + 



d9 2 



d9 



A 



m 2 2fi f (3 + 7 cos 9 
sin 2 fl ~ ~h? [ sin 2 9 



0(9) = 



(35) 



7 



(36) 



where m 2 and A = /(/ + 1) are the separation constants. The solution of Eq. ( |36l) is the 
well-known azimuthal angle solution. Eqs. ( 1341) and ( 1351) are the radial and the polar-angle 
equations. They have been examined separately in Refs. [3, In this paper, in order 
for the completeness of the paper, instead of just quoting their results and then discussing 
the effect of this new potential, we briefly show how to solve these equations by using the 
Nikiforov-Uvarov method {4] and then show the physical implications of this new non-central 
modified Kratzer potential. The radial part of the Schrodinger equation given by equation 
([34]) can be written as 



d 2 R lu 
+ 



dr 2 h 2 r 2 
This equation can be further arranged as 



(E - D)r 2 + 2Dar - Da 2 + 



Xh 2 
2/i 



R = 



(37) 



dr 2 

with the following abbreviations 

, 2n(E-D) 



sllR + l( e V-£r-K)f2 = 



(38) 



4fiDa ^ { D * 2 + f ) 

—^ir- ; K = tt. 



(39) 



h 2 ' * h 2 ' h 2 

It is now suitable for a NU solution. It is necessary to compare Eq. fl38l) with Eq. (JTJ) to find 
the solution of this equation. When these equations are compared, we obtain the following 
polynomials: 

f (r) = , a(r) = r , a{r) = e 2 r 2 — £r — k (40) 
If these polynomials are inserted into Eq. ([7j), we get n function as 



7T 



- ± ^f-Ae 2 r 2 + Ar(k + J) + 4k + 1 



(41) 



The expression in the square root must be the square of a polynomial according to the NU 
method. Thus, we can determine the constant k by using the condition that the discriminant 
of the square root is zero, that is, 



-f ± iey/l + 4k 



(42) 
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In view of that, new possible functions for each k are found as 

|±| [2re-iy/T+4Z\ , for k = -£ + is 

n = \ (43) 

\ ± | [2re + zVl + 4kJ , for k = -£- iey/1 + 4k 

To obtain the negative derivative of r = f + 27r, we select 



7r r 



- - ~ (W + zVTT4k") (44) 



and 



fc = -f - ie Vl + 4k (45) 
Using A = + tt' together with the values k and tt, r and A can be respectively obtained as 

r(r) = 1 + \/l + 4k - 2ire (46) 



A = -f - ze Vl + 4k - ze (47) 
Another definition of A n is given at Eq.flH]), 

A„ = 2me (48) 

and comparing this with Eq.( ]4"T{) . the exact energy eigenvalues of the radial part of the 
Schrodinger equation with the non-central modified Kratzer potential are derived as 

^ = D -^(l + 2, + ^l + 4(^ + ,( i + l)))" 2 (49) 

The energy eigenvalues calculated by this equation for the CO diatomic molecule are shown 
in the first column of Table [TT] for different n and I values, which are in agreement with the 
results of Ref. [8|. 

It is also possible to determine the radial eigenfunction of this potential by considering 
Eq.© and using Eq.gJ: 



0( r ) = r i/2(WI+3S) e -fcr (50) 

From Eq. (jSJ) and Eq. (jSJ), we get 

y (r) = ^L * ( r » p (r)) (51) 

yy ' p(r) dr n v Hy 11 y ' 

with p(r) = r Vi+4«; e -2i£r_ jf we g e ^ 2ier = C( r ) > \/l + 4k = the i? n ; radial wave function 
is 

iUO = ^C 1/2(1+ ^e- C/2 ^(C) (52) 



if the definition of U n i(r) = R n i(r)/r is used, the total unnormalized radial wavefunction is 

U nl (() = B*C xm - v) ^K{£) (53) 



where B„i is the normalization constant and it is evaluated as follows: 



B 



ni 



8fiDa 



3/2 



_(2n + C+l)(n + C)!. 



1/2 



(54) 



h 2 (2n + v + l) 

For the non-central modified Kratzer potential, the angle dependent Schrodinger equation 
is given in Eq. (1351) . We may also derive the eigenvalues and eigenfunctions of the polar 
angle part of the Schrodinger equation [7j similar to the solutions of the radial part. We can 
write Eq. ( 1351) by introducing a new variable |7|], x = cos 9, as follows: 

d 2 0(x) 2x d&(x) /A(l-x 2 )-m 2 -§ (P + jx) 



dx 2 



1 — x 2 dx 



+ 



(1-x 2 ) 2 



e(x) = o 



(55) 



To apply the Nikiforov-Uvarov method, Eq. fl55l) is compared with Eq. (TjQ) and the following 
polynomials are obtained: 



T 



-2x, a = 1 — x 2 , a 



-Xx 



> 2/^7 / 2 2fif3 

-W x+ [ x - m ~1F 



(56) 



The function 7r is obtained by putting the above expression in Eq. ([7]), 



k = ± \! x 2 (X — k) + -jrr- 1 ' -• ( A • - /;/- j-J- — k) 



(57) 



Thus, the polynomial of tt is found in four possible values: 



7T = ± < 



X 



X 



,2 , 3a3. 



,2 , 2ufl 



+ 



+ 



2 2 + 2i^- 



,2 i 2ufl 



, /or ft = 5 s yU 



(58) 



/or fc 



2\-m 



where u = 
derivative, 



m 2 _|_ _ (^I) 2 . For the polynomial of r = r + 2n which has a negative 



t = -21 



2 i 2h/3 
77T + -p 1 — U 



2x 1 + 



2 i 2u/3 , 



2 I V 2 

Using Eq. ([8]) and Eq. (Q, following expressions for the A are obtained respectively 

1 



(59) 



2A - (m 2 + 

A = - u — 

2 2 



m 2 + 2 $+U 



(60) 
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A n = 2n 1 + 



, m 2 + ^ + u 



+ n(n — 1) 



(61) 



If Eq. f )60|) and Eq. (16TT) are equated and the definition of A = 1(1 + 1) is used, the I values 
are obtained as 



+ n 



(62) 



If these I values are inserted into the eigenvalues of the radial part of the Schrodinger equation 
with the non-central Kratzer potential given by Eq. (1491) . the energy eigenvalues are found 
as follows: 



E„„ = D 



TP 



1 + In + 



( ( 

2/iDa 2 

fi2 



\ 


/ 




\ 


/ 


V 



(i.'i) 



We can also obtain the wave function of polar angle part of the Schrodinger equation [?| , 
using a and tt, 

4>{x) = (1 - x) B+c/2 (1 + xf~ c/2 (64) 



p{x) = (l-x*) B (\ + X 



X 



-c 



y n {x) = B n {\- x)- {B+C) (1 + x)- (B - C) 



dx r 



[ l + x f+^C {l _ x) n + B+C 



(65) 
(66) 



where B 



+ ^ + and C 



-2 i W 



The polynomial solution of y n is ex- 



2 V 2 

pressed in terms of Jacobi polynomials which are one of the ortogonal polynomials, giving 



P, 



(B+C, B-C) 



x). The corresponding wave functions are found to be 



0„(x) =N n (l- xY B+c)/2 (1 + xf B ~ c)l2 Pi B+c > B ~ c \x) 
where N„ is the normalization constant. 



(67) 



V. ENERGY EIGENVALUES USING AIM 

In this section, we show how to solve the Schrodinger equation for a particle in the pres- 
ence of non-central modified Kratzer potential by using AIM. In the spherical coordinates, 



11 



the equation is given by 

1 d ( 2 d\ 1 d ( . n d\ 1 d 2 2/i 



dr ^Yr ) + ^9 89 {^dO ) + r^TeW + f ^ ~ ^ ^' ^ = ° 

(68) 

If one assigns the corresponding spherical total wave function as ip(r, 9, (p)—-R(r)Y(9, (p), 
then by selecting Y(9, (p)= gin l 1/2e H(9)^(ip), the wave equation (IB"8"j) for a general non-central 
potential is separated into variables and the following equations are obtained: 

d 2 R{r) (In Kl + 1) 



dr , +(|(^-^))-^)i?(r)=0 (69) 

where m 2 and /(/ + 1) are separation constants. For bound states, we have the boundary 
conditions R(0)—0 and i?(oo)=0 in Eq. ([6"9l . H(0) and H(ir) are infinite in Eq. flTUj) and 
§((fi)=§((p + 2tt) in Eq. fl7ip . If we specialize to the case where V((p)=0, the normalized 
solution of the Eq. (ITTj) that provides the boundary condition is 

$ m (v>) = ^=e m \ (m = 0,±l,±2,...) (72) 

V 27T 



In this section we consider solutions of radial and angle- dependent parts of Schrodinger equa- 
tion for the non-central modified Kratzer potential within the framework of the asymptotic 
iteration method. The radial Schrodinger equation for the non-central modified Kratzer 
potential can be written as, 

« 2 «M + ( % ( E _ D ( , _^ + *)) _ >1±11) R(r) = o (73) 



2 \ f : jj^ 1 \ \ f ^"2 / / ^2 

If the following abbreviations are used: 

The radial Schrodinger equation takes the following form which is convenient in order to 
apply AIM. 



d 2 R{x) 
dx 2 



\ X X 2 J 



12 



d 2 R(x) + (_ e2 + 2a 2 _ r(r-l) \ 



<ix 2 



2 "" * /,'(./•) = 



x 



x z 



where 



T 



r-1) =a 2 + Z(Z + l) 



(76) 



(77) 



In order to solve this equation with AIM for I ^ 0, we should transform this equation to the 
form of Eg. (11 II) . Therefore, the reasonable physical wave function we propose is as follows 



R{x) = x T e- ex F(x) 



(78) 



If we insert this wave function into the Eg. (ITS"]) , we have the second-order homogeneous 
linear differential eguations in the following form 



d 2 F{x) 
dx 2 



x) dx 



x 



F{x) 



(79) 



which is now amenable to an AIM solution. To apply the AIM, it is reguired to compare 
Eg. (1791 with Eg. (122]) . Subseguently, by using Eg. ([TT]) . the values of X n (x) and s n (x) are 
computed as follows 



Ao(a;) 
s (x) 

Ai(x) 

Si(x) 

A 2 (x) 
. . . etc 



2( £ -I 

X 

2 



x 
2t 



(er - a 2 ) 



2t 



— + - (er - a 2 ) + 2e 

X A X \ X 



/2e It + 1 
2(--« 2 )(-- — 



4 



a — er 



2e - 



2r 



X" V X 

4er - 2a 2 2er - 2a 2 



6r 4(er — a 2 ) 

X z X 



X 



6r 2(e(r-l 
h — — 

2 ~ 



x 



x 



2e 



2e 2r Y 
x 2 J 



2r 



(80) 
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If we use the termination condition of the AIM given in Eq. (fl9l) . energy eigenvalues are 
obtained as follows 



£o _ £1 
Ao Ai 

£l _ £2 

Ai A2 

S2 _ S3 

A2 A3 



a 

=> e = — 

r 



a 



62 



r+ 1 

,,2 



a 



t + 2 
. . . etc 



(81) 
(82) 
(83) 



which can be generalized as 



a 



(84) 



t + n 

If one inserts the values of e, r, and a into equation ( 1841) . the ro- vibrational energy spectrum 
of the Schrodinger equation with the non-central modified Kratzer potential becomes 

-2 



E. 



Til 



D 



h 2 



l + 2n 



1 + 4 |'^ + /(i + i : 



(85) 



where n and / are the vibrational and the rotational quantum numbers, respectively. 

As indicated in Section IHIt we can construct the corresponding eigenfunctions by using 
Eq.pPI). Eq.(TT9l) provides Eq.j22J for N=-l, b=0, a=e and m=r-l. 

Equation ( 13T1) can be written using limit relation as b — > in the following form 



\im 2 Fi(—n, lib + a; c; z6) = i-FU— n; c: z) 



(86) 



Vn{X) 



■l) n C 2 (N + 2r(a) nl F 1 (-n,a; 



2a 



-x 



N+2\ 



N + 2 

Directly solutions F n (x) for the eigenvalue problem Eq. fl85|) can be obtained from Eq.([8 
with the substitution N=— 1 and a=e. 



(87) 



F n {x) = {-l) n C 2 (a) n iFx (-71, a; 2ex) 



(88) 



where <j=2t, and (cr) n 



r(cr+n) 



r ^p. Consequently, the radial eigenfunction for Schrodinger 
equation with modified Kratzer potential is 



R( x ) = x T e~ €X (-1)™ C 2 (a) n t Fi (-n, 2r; 2ex) 



R{r) 



r e~^ (-l) n C 2 (a) n x Fi (-n,2r; 2e^) 



(89) 
(90) 
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where C2 is normalization constant. 

We now probe solutions of the angle-dependent Schrodinger equation with the non-central 
modified Kratzer potential using the same approach. The non-relativistic angular motion of 
a diatomic molecule of mass \i is described by the following equation, 



d 2 H(9) (2n 



d6 2 



m 



1 1 cos 2 6 



t9 Vi ~ d) + sm 2 6 2 Asm 2 9 



-1(1 + 1) H(6)=0 



(91) 



If we define as l=l'-\ and 1(1 + 1)=/ /2 -^ , Eq. (l9~Tj) can be written as 



d 2 H(9) 
dd 2 



2fip/h 2 + m 2 - 1/4 2yrfcos9/ti 



+ 



sin 2 6 sin 2 6 

This equation can be further arranged as 

d 2 H(6) \k 2 + t] 2 -1/A 2kt)cos9 



H(6) = -l' 2 H(6) 



(92) 



with 



d6 2 



K 



sin 2 d 



+ 



sin 2 ^ 



H(6) = -l' 2 H(6) 



2tf 

h 2 



+ m + \ (m 2 + 



h 2 ' 



2^7 
h 2 



W M 2 



(m 2 + ~0-) 2 



2^7 

h 2 



(93) 



(94) 



(95) 



Equation ( J93l) can be written in the following form by introducing a new variable of the 
form x = cos8, 

I' 2 



d 2 H(x) x dH(x) 



dx 2 



1 — x 2 dx 



K + TTj + 2kT]X 



(l-x 2 ) 2 



H(x] 



l-x 2 



;H(x) 



(96) 



Let the angular wave function be factorized as: 



H(x) = (1 - x)^ +1 / 4 \l + xf-^ +1 '^ f(x) 



(97) 



Equation (I9"6"|) reduces to the second-order homogeneous linear differential equation in the 
following form 

d 2 f(x) 



dx 2 



2 [ V + (k + 1)x } df(x) + ( ( K + l/2) 2 -l' 2 ^ f( r) 



l-x 2 



dx 



l-x 2 



(98) 



which is convenient to a AIM solution. In order to find the solution of this equation, it is 
necessary to compare Eq. (|98l with Eq. tflT!) . By means of Eq. (TlTI) . the values of X n (x) and 



15 



s n (x) are obtained as follows 
2r/ + (2k + 2)x 



Ao(ar) 
s Q (x) 
Ai(x) 

si(x) 

A 2 (x) 



etc 



1-x 2 

\2 j/2 



(k + 1/2)"-/' 
1 — X 2 

(k + 1/2) 2 - /' 2 + 2/t + 2 (2?7 + (2k + 2)x) (2rj + (2k + 4)x) 



(1-x 2 ) 



2^ 



((k+1/2)^-/ /2 )(2t ? + (2k + 4)x) 
(1-x 2 ) 2 

2(t/ + (k + 1)x)((k + 1/2) 2 - I' 2 + 6k + 8) 
(1-x 2 ) 2 

2(r/ + (k + 1)x)((k + 1/2) 2 - I' 2 + 8x 2 + (2r/ + (2k + 2)x)(2r/ + (2k + 8)x)) 



2A3 



(1-X 

( (k + 1/2) 2 - / /2 )((k + 1/2) 2 - V 2 + 2(2k + 2)) 
(1-x 2 ) 2 

((k + 1/2) 2 - Z ,2 )(8x 2 + 6(2r/ + (2k + 2)x)x + (2r/ + (2k + 2)x) 2 ) 
+ (1-x 2 ) 3 



(99) 



Combining these results with the condition given by Eq. (|T9|) yields 

f =► I' 2 ' 1 
Ai 



so 
Ao 

£i 
Ai 

A 2 
. . etc 



S2 

A 2 

S3 

A 3 



K + 



/2 



7/2 



K + 



K + 



(100) 
(101) 
(102) 



When the above expressions are generalized, the eigenvalues turn out as 

2 



l2 



Inserting k and V in Eq. (11031) 



k + N + - 
2 



iV = 0,1,2, .. 



(103) 



\ 



^ + m 2 + /( m 2 + ^ 



'-M/3 ^ 2 _ ( Vl)2 



(104) 



If / value obtained by Eq. (11041) is inserted into energy spectrum of radial part of the 
Schrodinger equation given by Eq. (l85l) . we find energy spectrum for a diatomic molecule 
system in the presence of non-central modified Kratzer potential as following, 
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This eigenvalue equation obtained by using the asymptotic iteration method is the same 
as the EH1 As similar to radial wavefunction, angular wavefunction can be found by using 
Eq.( l30i) . The advantage of the asymptotic iteration method is that it gives the eigenvalues 
directly by transforming the second-order differential equation into a form of y" =Xo{r)y' + 
s (r)y. The wave functions are easily constructed by iterating the values of s (r) and A (r). 
The method presented in this study is a systematic one and puts no constraint on the 
potential parameter values involved. 

VI. INTERPRETATION OF THE RESULTS 

In comparison with the eigenvalue equation (H9l) of the modified Kratzer potential, the 
eigenvalue equations given by equations ([63] and I105f) have the correction due to the angle- 
dependent (non-central) part of the modified Kratzer potential, given by equation ( |32l ). 
This correction has two dependencies: First one is due to the values of the constants (3 
and 7 in the non-central part of the potential. Different values of the /3 and 7 change the 
shape and the depth of the potential which affects the values of the energy eigenvalues. The 
second one is because of the orbital angular momentum quantum number I given by equation 
( |62i) . It is known that the orbital angular momentum quantum number / is a good quantum 
number and it is a constant of motion for the central potentials such as the modified Kratzer 
potential. On the other hand, for the non-central modified Kratzer potential, it is no longer 
a constant of motion and its value is determined in terms of the constants f3 and 7 as well 
as n and m quantum numbers as given in equation (|62|) . This is the effect of the non-central 
part of the potential which modifies the I values. For /3=7=0, the non-central modified 
Kratzer potential reduces to the modified Kratzer and they should give the same energy 
eigenvalues for the same I values. For the modified Kratzer potential, / is a constant of the 
motion and takes the values Z=0,l,2..., but for the non-central modified Kratzer potential, it 
is determined from equation (I62I) . Therefore, they give the same energy eigenvalues as long 
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as the following condition is satisfied: the / value of the modified Kratzer potential is equal 
to (n + m) of the non-central modified Kratzer potential for [3—j—O. By comparing the 
eigenvalue equation (E nm ) for the non-central Kratzer potential given by equations (16Uland 
11051) with the eigenvalue equation (E n i) for the modified Kratzer potential given by equation 
( 149 p . we note that the effect of the non-central part is small for the small values of f3 and 
7 constants, but the difference becomes apparent when they have large values. In order to 
show this, the difference between the eigenvalues of the modified Kratzer and non-central 
modified Kratzer potentials is demonstrated in comparison with each other in Table [III for 
different n, I and m quantum numbers. 

VII. CONCLUSION 

In this study, we have calculated the exact bound-state energy eigenvalues and the corre- 
sponding eigenfunctions of the new exactly solvable non-central modified Kratzer potential 
by using two different methods. Both NU and AIM methods generate the same results. 
In comparison with the NU method, AIM puts no constraint on the potential parameter 
values involved. Therefore, Aim is more systematic than the NU method in solving such 
second-order differential equations. 

In our study, we have determined the effect of the non-central term on the bound-state 
energy eigenvalues explicitly. The energy eigenvalue equations for the modified and non- 
central modified Kratzer potentials are given by equation (j4"9~l) and equations (1631 and ["1051) . 
respectively. In comparison with the eigenvalues of the modified Kratzer potential, the effect 
due to the angle-dependent part of the non-central potential on the energy eigenvalues is 
shown for the CO molecule in Table [TTJ The correction is very small for the small values 
of 13 and 7 in the angle dependent part of the potential. However, the difference becomes 
apparent for the values of /3=j=l and over. In order to examine this effect, we have plotted 
the shape of the modified Kratzer and non-central modified Kratzer potential in Figure [1] 
for different values of f3 and 7. We perceive from this figure that the sum of the modified 
Kratzer and centrifugal potentials with 1=50 gives almost the same shape for the non-central 
modified Kratzer potential with /5=7=1. That is, the angle-dependent part behaves like the 
centrifugal barrier and as it is clearly seen in Figure [lj it reduces the depth of the potential 
pocket, which effects the bound-state energy eigenvalues. 
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FIG. 1: Comparison of the modified Kratzer and non-central modified Kratzer potentials {6 = 
for different (3 and 7 values for the CO diatomic molecule. 
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30°) 



Molecule D e (eV) a (in A ) fx (in amu) 



CO 10.84514471 1.1282 6.860586 

TABLE I: Reduced mass and spectroscopical properties of the CO diatomic molecule in the ground 
state 

TABLE II: Comparison of the eigenvalues of the modified Kratzer (E n i^) and non-central mod- 
ified Kratzer (-Enm(Arc-A')) potentials for different n, I and m values with (5 = 7=0.0, 0.1, 1.0 and 
5.0 values for the CO diatomic molecule, calculated by using equations (f4U|) . ([63]) and (|105p . 
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